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COULOMB  FRICTION,  PLASTICITY.  AND  LIMIT  LOADS1 

by  D.  C*  Drucker^ 


Abstract 

Additional  attention  is  given  to  the  somewhat  subtle 
but  extremely  important,  difference  between  Coulomb  friction  and 
the  apparently  corresponding  resistance  to  plastic  deformation. 
It  is  shown  that  the  limit  theorems  previously  proven  for 
assemblages  of  perfectly  plastic  bodies  do  not  always  apply 
when  there  is  finite  sliding  friction.  Theorems  are  developed 
which  relate  the  limit  loads  with  finite  Coulomb  friGtion  to  the 
extreme  cases  of  zero  friction  and  of  complete  attachment,  and 
also  to  the  case  where  the  frictional  interfaces  are  n cemented51 
together  with  a cohcsionless  soil. 


The  results  presented  in  this  paper  were  obtained  in  the  course 
of  research  sponsored  by  the  Office  of  Naval  Research  under 
Contract  N7onr-35801  with  Brown  University, 


Professor  of  Engineering,  Brown  University,  Providence  12,  R.  I. 
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Sliding  Friction  vs.  Plastic  Resistance » 

3 

It  was  pointed  out  previously  that  the  similarity 
between  the  stress-strain  relation  for  a rigid-  or  elastic- 
perfectly  plastic  material  and  the  force-displacement  relation 
for  the  Coulomb  friction  case,  Fig.  1,  may  be  misleading.  The 
essential  feature  may  best  be  stated  in  mechanical-thermodynamic 
terms  for  an  assemblage  of  elastic-plastic  bodies  in  equilibrium, 
Fig.  2,  if  either  the  attachment  is  complete  or  if  all  coef- 
ficients of  friction  are  zero,  no  work  can  be  extracted  from 
the  bodies  and  any  equilibrium  system  of  forces  acting  upon  them. 
In  other  words,  there  is  no  way  in  which  the  bodies  plus  the 
forces  can  act  as  an  engine  in  this  thermodynamic  sense.  They 
can,  however,  if  a coefficient  of  friction  is  finite^. 

This  point  of  view  is  made  clearer  and  less  abstract 
by  a comparison  between  Coulomb  sliding  friction  and  shearing 
on  a plane  through  a Coulomb  cohesionless  soil  obeying  the 
generalized  plastic  potential  laws,  Fig.  3*  The  relation  between 
the  normal  vortical  force  N and  frictional  horizontal  force 
F In  each  case  may  be  written  as  F ==  p,  N.  The  displacement 
picture,  however,  is  fundamentally  different^.  The  block  In 
the  friction  case  slides  in  the  direction  of  the  force  F.  It 
also  moves  up  in  the  ease  of  shearing  of  soil  because  of  the 
volume  expansion  which  accompanies  shear  if  the  generalized  plas- 
tic  potential  relations  are  followed  . The  displacement  or 

^D.C,  Drucker,  "Some  Implications  of  Work-Hardening  and  Ideal 
Plasticity" , Quarterly  of  Applied  Mathematics,  v.  7 (1950) 
pp.  4-11-418. 

4. 

D.C.  Drucker  and  W.  Pragcr,  "Soil  Mechanics  and  Plastic  Analysis 
or  Limit  Design",  Quarterly  of  Applied  Mathematics,  v.  1C  (1952) 
PP.  157-165. 
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velocity  vector  makes  an  angle  cp  with  the  horizontal  where 

-1 

9 = tan  M-  is  the  angle  between  the  straight  line  envelope  to 
the  limiting  Mohrs  circles  and  the  negative  d axis.  Work  is 
done,  therefore,  against  the  downward  force  N.  It  is  this 
negative  work  which  prevents  elastic-plastic  bodies  and  the 
system  of  forces  acting  on  them  from  acting  as  an  engine  while 
apparently  similar  frictional  systems  can. 

Limit  Theorems. 

5 6 

The  two  major  limit  theorems  * in  a crude  sense  state 
that  an  assemblage  of  elastic-perfectly  plastic  bodies,  with 
zero  friction  or  complete  attachment  at  each  interface,  will  on 
the  one  hand  do  the  best  they  can  to  distribute  stress  to  avoid 
collapse  and  on  the  other  will  recognize  defeat  if  any  kinematic 
collapse  mode  exists.  This  anthropomorphic  approach  is  refined, 
when  geometry  change  is  negligible,  in  the  actual  statement  of 
the  theorems  proved  as: 

1.  Collapse  will  not  occur  if  any  state  of  stress  can  be 
found  which  satisfies  the  equations  of  equilibrium  and  the  boundary 
conditions  on  stress  and  which  is  "below  yield"  at  each  point. 

2.  Collapse  must  occur  if  for  any  compatible  flow  pattern, 
considered  as  plastic  only,  the  rate  at  which  the  external  forces 
do  work  on  the  bodies  equals  or  exceeds  the  rate  of  internal 
dissipation. 

3.  Collapse  takes  place  at  constant  stress  so  that  strain 
rates  are  purely  plastic. 

^D9C.  Drucker,  H.J,  Greenberg,  M,  Prager.  "The  Safety  Factor  of 
an  Elastic-Plastic  Body  in  Plane  Strain",  Journal  of  Applied 
Mechanics,  v.  18,  Trans.  ASME,  v.  75  (1951)  pp.  371-378. 

D.C.  Drucker,  H.J.  Greenberg,  V/,  Prager,  "Extended  Limit  Design 
Theorems  for  Continuous  Media",  Quarterly  of  Applied  Mathematics, 
v.  9 (1952)  pp  381-389, 
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Nov;  suppose  that  the  coefficient  of  friction  is  finite 
and  non-zero  on  at  least  one  of  the  interfaces  S of  Fig.  2. 

Is  the  system  still  intelligent  enough  to  distribute  the  stress 
to  avoid  collapse,  and  arc  the  theorems  still  valid?  Trouble 
immediately  appears  in  the  upper  bound  Theorem  2.  The  rate  of 
internal  dissipation  cannot  even  be  calculated  in  all  cases  because 
frictional  dissipation  is  not  determined  uniquely  by  the  flow 
pattern.  It  depends  not  alone  upon  relative  displacement  rates 
but  also  on  the  normal  pressure  on  the  frictional  interface5  a 
quantity  which  will  often  not  bo  known.  This  type  of  difficulty 
does  not  appear  in  the  lower  bound  Theorem  1.  It  might  seem 
plausible  to  assume  that  the  theorem  is  valid  with  the  additional 
requirement  that  the  state  of  stress  not  violate  the  friction 
condition  at  the  interface  as  well  as  staying  below  yield.  A 
simple  example  shows  this  hopeful  intuitive  approach  to  be  in 
error. 

Illustrative  Examples. 

Fig.  4a  represents  two  rigid  blocks,  the  heavy  horizontal 
one  is  balanced  on  a rough  ledge  while  the  light  small  one  is  on 
an  inclined  plane  with  friction  angle  tp  and  inclination  at.  The 
contact  surface  between  the  blocks  is  at  a slightly  flatter  angle 
P than  the  plane  and  will  be  considered  frictionless.  There  need 
bo  no  force  between  the  blocks  and  yet  a large  force  N eould  be 
carried  between  the  blocks  without  violating  equilibrium, Fig.  4b. 

If  cp  > a-fi  , such  a force  is  stabilizing.  Clearly,  however, 
if  a > cp  the  block  will  slide  dow;i  the  plane  and  no  stabilizing 
force  will  be  developed.  On  the  other  hand,  if  the  sliding  was 
really  plastic  shearing  as  for  a soil,  then  the  incipient  velocity 
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vector  V would  point  away  from  the  inclined  plane  at  an  angle 
<p  as  in  Fig.  4c,  The  large  block  would  have  to  lift  up  and  the 
force  N would  indeed  be  mobilized  to  prevent  sliding  down  the 
plane.  A modified  Theorem  1 is  thus  seen  to  be  improper  for 
the  friction  case  but  is  again  verified  for  the  ’•plasticity" 
problem. 

The  problem  of  Fig.  4 might  seem  tricky  and  exceptional 
because  the  little  block  is  not  really  confined  by  the  large  one 
and  because  the  bodies  are  rigid,  A more  elaborate  example  is 
probably  needed  to  demonstrate  the  point  properly.  The  punch  or 
equally  well  the  uniform  pressure  example  of  Fig,  5 has  special 
dimensions  to  avoid  calculations  but  is  otherwise  a good  repre- 
sentative plane  strain  problem.  The  interface  UMNW  between 
the  upper  arch-shaped  block  and  the  lower  support  is  supposed  to 
have  an  extremely  large  but  finite  friction  angle.  The  yield 
stress  in  shear  will  be  called  k and  either  the  Mises  or  Tresca 
yield  condition  is  employed. 

Fig.  shows  a discontinuous  stress  solution  passing 

through  the  interface  as  though  the  two  bodies  were  one.  This 

o 7 

is  the  often  used  30  wedge  solution  slightly  modified  , In  the 

equilateral  triangle  region  ABC  the  principal  stresses  are  com- 
pressive and  of  magnitude  3k  in  the  vertical  direction  and  k 
horizontally.  In  the  sloping  region  BCDE  the  stress  is  2k 
parallel  to  BE  and  zero  in  the  perpendicular  direction.  Region 
DHE  has  principal  stresses  +k  and  -k.  The  shearing  stress, 
if  any,  and  the  normal  traction  on  the  pianos  of  discontinuity 

?P,  G.  Hodge  and.  W,  Prager,  "Theory  of  Perfectly  Plastic  Solids", 
John  Wiley  & Sons,  1951,  Section  26, 
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BC,  DE,  HE,  HD  are  continuous  across  the  planes. 

Both  compressive  and  shearing  tractions  cs  and  t in 
magnitude  which  are  easily  calculated,  will  act  on  the  interface 
UMNW.  The  coefficient  of  friction  p is  assumed  large  enough 
so  that  T<p  c , Therefore,  the  stress  field  of  Fig.  5 a multiplied 
by  a number  ever  so  slightly  less  than  unity  will  be  "below  yield" 
everywhere  and  does  not  violate  the  friction  condition.  The  total 
force  P per  unit  dimension  perpendicular  to  the  paper  can  be 
made  as  close  to  3kb  as  desired. 

The  question  of  whether  or  not  the  bodies  are  able  to 
adjust  themselves  in  this  static  manner  is  answerdd  by  the  possible 
discontinuous  rigid-block  kinematic  collapse  pattern,  Fig.  5b. 

Block  ATB  is  taken  to  more  vertically  downward  with  velocity 
V requiring  ATMU  and  BTNW  to  move  at  45°  with  velocity  V/V? 
so  that  separation  occurs  on  MU  and  NW.  Block  MTN  is 
stationary.  The  velocity  discontinuity  on  the  surfaces  of  sliding, 
3TM  and  ATN  each  of  length  b/2  , is  V/72.  The  rate  of  energy 
dissipation  per  unit  perpendicular  dimension  is  just 
2Ck  • oi f'l  * V//2)  = 2kbV  because  in  this  permissible  virtual 
displacement  pattern  there  is  separation  along  the  friction  surfaces 
MU  and  MW  and  they  contribute  nothing.  Equating  the  rate  of 
work  of  the  external  forces  PV  to  2kbV  gives  the  upper  bound 
and  in  this  case  correct  result  for  Fig.  5,  P = 2kb,  and  not 
b £ ■ ■ • 

The  intuitive  limit  theorem  on  the  inherent  intelligence 
of  the  material  is  thus  seen  to  be  incorrect  for  finite  and  non- 
zero friction.  Again,  substitution  of  an  elastic-plastic  "soil" 

for  the  Coulomb  friction  interface  fixes  matters.  Surfaces  MU 

If 

and  NW  do  increase  the  energy  dissipation  and  the  value  3kb 
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is  a proper  lower  bound  if  the  angle  <p  is  large  enough. 

It  might  be  thought  that  the  separation  counter-examples 
given  in  this  section  are  unfair.  However,  it  is  precisely 
the  freedom  in  the  friction  case  to  slide  with  an  arbitrary 
normal  stress  on  the  surface  of  sliding  and  therefore  an 
arbitrary  dissipation  which  prevents  the  application  of  the 
11  intuitive  theorems”.  Sliding  at  zero  normal  stress  and  actual 
separation  are  essentially  the  samo.  An  angle  of  45°  instead  of 
44°  in  Fig.  5 would  not  really  change  the  result. 

Friction  Theorems. 

The  question  now  arises  as  to  what  can  be  said  about 
limit  loads  for  assemblages  of  bodies  with  frictional  interfaces. 
Two  theorems  seem  intuitively  obviousi 

A,  Any  set  of  loads  which  produces  collapse  for  the  condi- 
tion of  no  relative  motion  at  the  interfaces  will  produce 
collapse  for  the  case  of  finite  friction.  No  relative 
motion  is  a more  inclusive  term  than  infinite  friction 
because  separation  is  not  permitted, 

B,  Any  set  of  loads  which  will  not  cause  collapse  when  all 
coefficients  of  friction  arc  zero  will  not  produce 
collapse  with  any  values  of  the  coefficients. 

Although  at  this  stage  intuition  may  not  seem  completely 
reliable,  the  theorems  are  in  fact  true.  Their  proof  follows  in 
part  the  technique  employed  in  developing  the  limit  theorems  for 
plastic  bodies^1.  The  most  important  tool  is  the  theorem  of 
virtual  work 

J'ATiVA+  iTFl5idV  =Iv°lfi3dV  + Is  TiA“idS 
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in  which  the  summation  convention  is  employed  . There  is  no 
necessary  tie  between  the  equilibrium  system  of  surface  tractions 
Tj_,  stresses  °ij  > and  body  forces  F^,  and  the  compatible 
system  of  velocities  u^,  strain  rates  and  velocity  dis- 

T 

continuities  Au^.  The  surface  area  A of  the  assemblage  of 

bodies  of  volume  V does  not  include  the  frictional  interfaces 

S.  For  convenience  of  description,  the  velocity  will  be  assumed 

continuous  except  across  S.  The  term  for  velocity 

discontinuities  in  the  material  is  omitted  but  this  does  not 

5 6 

actually  restrict  the  generality  of  the  result*  * 

Use  will  also  be  made  of  the  properties  of  the  yield 

p 

surface  f(c4J)  = k . Fig.  6,  The  rate  of  dissipation  of  energy 
i j 


• p 

per  unit  volume,  D , 


is  uniquely  determined  by  the 


plastic  strain  rate  e^jj.  T*10  ,,vector't  e.^,  Fig.  6,  is  normal 


ij 


to  the  yield  surface  at  a smooth  point  or  lies  between  the  normal 
to  the  surface  at  adjacent  points  to  a vertex  or  corner.  Further- 
more, the  surface  is  convex  so  that  the  dot  product  of  the  plastic 

i 


strain  rate  ’’vector"  with  any  stress  vector  to  a point 

f 2 

the  surface,  f(o^j)  < k , cannot  be  as  large  as  D, 


ij 


inside 


- i • p • p 

Vu  < D - Vu 


/ O \ 
\ / 


,N  N 


To  prove  Theorem  A,  assume  that  a set  of  loads  Tj.,  F^ 
produces  collapse  when  all  relative  motion  is  prevented  at  the 
interfaces  S (bodies  welded  together).  This  type  of  composite 
body  can  be  analyzed  by  the  established  limit  theorems  which  state 


mN*N, . I _H*N.  x 

T-iUjdA  + F^Ujdv  >. 
A Jv 


n N*  N 

°ij£ijdv 


(35 
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Reference  7,  Chapter  8,  As  an  example, 

*T  *-r  *T 

= Vx  + aycy  + dzsz  + TxyY  xy 


yz'  yz 


* T 

+ TzxYzx 


I 


1 
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' N 
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is  purely  plastic, 

• N 


e ^ (elastic)  = 0 


(4) 


For  Theorem  A to  be  false,  there  must  exist  some  equilibrium 

A 


state  of  stress  cJ 


13 


such  that  f(c^j)  < k^.  Using  the  virtual 
N * n 


work  Equation  (1),  with  Uj_,  e ^ as  the  compatible  state, 


r 


T?u?dA  + I F?u?dv  = 

JA  11  Jv  1 i 


*N  _ 


A*  N 

6 ije  13dv  + 


TjAuNdS 


(5) 


However,  Au  = 0 because  the  velocity  solution  represents  no 
relative  motion  at  S and,  from  Inequality  (2),  < di^ei^* 

Equation  (5)  is,  therefore,  in  contradiction  to  Inequality  (3) 
and  Theorem  A can  not  be  false. 

The  proof  of  Theorem  B starts  from  the  given  set  of 
loads  T^,  and  stresses  which  do  not  cause  collapse 

with  p.  = 0.  Next  suppose  that  collapse  can  occur  for  \i>  0 
with  a collapse  field  uf,  e i j j Au^  and  stress  field 
Virtual  work  Equation  (1)  becomes 

J.  T^dA  + J v F^dv  = L + L TKas  (6) 


and  the  collapse  condition  is 

f T?u?dA  + | F?u?dv  > 

JA  Jv  1 1 ~ J 


„ B* 

0 j * £ -t  1 dv  + 


ij  13 


T^AufdS 


(7) 


Conditions  (6)  and  (7)  are  incompatible  and  Theorem  B cannot  be 

B*  B 0*  B B *B 

false  because  > d^e^,  T^Au^  is  a frictional  dissipa- 

0 *B 

tion  and,  therefore,  zero  or  positive,  and  T^Au^  is  zero  or 


negative.  The  last  statement  follows  from  the  fact  that  T 


0 


• T> 

is  normal  to  S at  each  point  and  is  compressive  while  Au£  is 
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a relative  tangential  displacement,  and  possibly  a separation  as 
well,  but  not  an  overlap. 

The  friction  theorems  A and  B have  thus  been  proved  for 

all  stable  convex  yield  functions.  They  occasionally  enable  the 

limit  load  to  be  computed  quite  precisely  for  finite  non-zero 

friction.  The  well  known  two  dimensional  punch  problem  for  a 

Prandtl-Reuss  or  a Mises  material  provides  such  an  example.  Two 

solutions  are  available  for  upper  bound  computations.  One  by 
9 

Prandtl  , Fig.  7a,  contains  a rigid  region  which  acts  as  an  ex- 
tension of  the  punch;  there  is  no  relative  motion  between  the 
punch  and  the  contact  area.  The  other  by  Hill‘d,  Fig.  7b,  assumes 
Eero  friction  and  appreciable  slip  does  take  place.  Both  solu- 
tions give  the  same  answer  for  the  average  pressure,  p = (2+it)k 
■where  k is  the  yield  stress  in  shear.  A lower  hound  solution 
of  5k  has  been  obtained”^'.  Therefore,  the  limit  pressure  is 
between  5k  and  5.14k  for  all  possible  values  of  the  coefficient 
of  friction. 

It  will  often  be  found,  however,  that  Theorems  A and 
B do  not  provide  very  close  bounds.  The  concept  of  the  plastic 
cohesionless  soil  interface,  already  discussed  in  considerable 
detail,  may  then  be  of  furthci  help  in  the  following  theorem, 

9 

L,  Prandtl,  "Uaber  die  Haerte  plastischer  Koerper" , Nachrichten 
von  der  Koeniglichen  Gesellschaft  der  Vissenschaf ten  zu 
Goettingen,  iiathematisch-physikalische  Klasse  (1920)  pp.  74-85. 

^R*  Hill,  “The  Plastic  Yielding  of  Notched  Bars  under  Tension", 
Quarterly  Journal  of  Mechanics  and  Applied  Mathematics,  v.  2 
(19^9)  PP.  40-52. 

11 

R.  T.  Shield  and  D.  G.  Drucker,  "The  Application  of  Limit  Analy- 
sis to  Punch  Indentation  Problems",  to  appear  in  Journal  of 
Applied  Mechanics, 
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C.  Any  set  of  loads  which  w 111  not  cause  collapse  of  an 
assemblage  of  bodies  with  frictional  interfaces,  will 
not  produce  collapse  when  the  interfaces  are  "cemented" 
together  with  a cohesionless  soil  of  friction  angle  <p  = 
arc tan  ^ . 

The  proof  of  Theorem  C follows  essentially  from  the 

observation  that  the  state  of  stress  in  the  friction  case  satisfies 

2 

the  conditions  < k and  |x|  _^p,|c|at  impending  collapse. 

A safe  state  of  stress  exists,  therefore,  for  the  soil  case  when 

collapse  does  not  occur  in  the  friction  problem. 


Conclusions 

The  limit  load  for  an  assemblage  of  bodies  with 
frictional  interfaces  is  bounded  below  by  the  limit  load  for  the 
same  bodies  with  zero  friction  on  the  interfaces.  It  is  bounded 
above  by  the  limit  load  for  no  relative  motion  at  the  interfaces 
and  also  by  the  limit  load  for  the  same  assemblage  cemented  at 
the  interfaces  by  a cohesionloss  soil.  Limit  theorems  applicable 
to  these  bounding  problems  do  not  apply  generally. 
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FiG.  2.  AN  ASSEMBLAGE  OF  BODIES 
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FIG.  4.  A COUNTER-EXAMPLE- RIGID  BODIES 


FIG.  5.  TWO  DIMENSIONAL  PROBLEM  WITH  ELASTIC- PLASTIC 
BODIES  (a)AN  EQUILIBRIUM  SOLUTION  lb)  A VELOCITY  (KlNEMATiCJ 

SOLUTION  INVOLVING  SEPARATION 
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FIG.  6.  THE  YIELD  SURFACE 


(a)  (b) 


FIG.  7.  (a)  PRANOTL  AND  lb)  HILL  SOLUTIONS  TO  PUNCH  PROBLEM 


